Observation of negative absolute resistance in a Josephson junction 
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We experimentally demonstrate the occurrence of negative absolute resistance (NAR) up to about 
— Ifi in response to an externally applied dc current for a shunted Nb-Al/A10i,-Nb Josephson 
junction, exposed to a microwave current at frequencies in the GHz range. The realization (or not) 
of NAR depends crucially on the amplitude of the applied microwave current. Theoretically, the 
system is described by means of the resistively and capacitively shunted junction model in terms 
of a moderately damped, classical Brownian particle dynamics in a one-dimensional potential. We 
find excellent agreement of the experimental results with numerical simulations of the model. 

PACS numbers: 05.45.-a, 05.40.-a, 05.60.Cd, 74.50. +r 



When a static force is applied to a system consisting 
of mobile particles, these particles usually move in the 
direction of the force, i. e., they show positive mobility, 
which leads to, e. g., a positive conductance or resistance 
in electrical systems. Also well known is the fact that 
such a system can exhibit regions of negative differential 
mobility/resistance [l|, [3, S Ij, 0, Q- However, the ab- 
solute mobility/resistance usually remains positive. The 
opposite response, i. e., a motion against the static force 
is termed negative absolute mobility or negative absolute 
resistance (NAR) . This is clearly a quite counter intuitive 
effect which, at first glance, might seem even to be in con- 
flict with Newton's laws and thermodynamic principles 
Q. Yet, nonlinear systems being driven far from equi- 
librium can indeed exhibit not only a negative differen- 
tial resistance but also a NAR effect. Unambiguous and 
convincing experimental observations of NAR are still 
quite scarce, involving systems consisting of electrons in 
a sample of bulk GaAs Q, electrons in semiconductor 
heterostructures Q , electrons in low dimensional conduc- 
tors (To| . and charged Brownian particles in structured 
microfluidic devices pd| . Apart from the low dimensional 
conductors, the system was always driven out of equilib- 
rium by means of an ac driving force and then its re- 
sponse to an externally applied static perturbation was 
studied. On the theoretical side, a considerably larger 
literature is available, most notably on different types of 
semiconductors and semiconductor heterostructures 
In all those cases (except [ll|) NAR is based on purely 
quantum mechanical effects which cannot be transferred 
into the realm of classical physics. For classical systems, 
a first theoretical demonstration of the effect was pro- 
vided in the context of a spatially periodic and symmet- 
ric model system of interacting Brownian particles, sub- 



jected to multiplicative white noise While each of 

the different ingredients of the model is quite realistic 
in itself, their combined realization in an experimental 
system seems difficult. In particular, the main physical 
mechanism is based on collective effects of at least three 
interacting particles [isj . An entirely different mecha- 
nism was later on suggested theoretically for a realistic, 
classical model dynamics of a single Brownian particle in 
a suitably tailored, two-dimensional potential landscape 
in Ref . 1 1 411 a nd subsequently realized experimentally in 
Refs. [Ill, [131 . As a first application of NAR, the separa- 
tion of different particle species has been realized in Ref. 
[l6t . While the underlying basic physical mechanism still 
requires at least two spatial dimensions, very recently, 
NAR has been analyzed and predicted theoretically to 
occur also in the simplest possible case of a single Brow- 
nian particle dynamics in one dimension Ell • More 
precisely, two basically different physical mechanisms ca- 
pable of generating NAR in such systems have been un- 
raveled, namely a purely noise induced effect in Ref. [l7j 
and a transient chaos induced effect in Refs. [ll]. In both 
cases, an experimental realization by means of a Joseph- 
son junction subjected to suitable dc and ac currents has 
been proposed. In this Letter we show that a moderately 
damped Josephson junction being driven by microwaves 
indeed shows NAR of the type predicted in Refs. [ist . 
A first hint along these lines can be found in Fig. 13 of 
p^ . although without further explanation or discussion 
and no direct reference to the resistively and capacitively 
shunted junction model. 

To model the Josephson junction we use the resistively 
and capacitively shunted junction model [13, HH- It de- 
scribes the equation of motion for the difference 6 of the 
phases of the superconducting order parameter in the two 
electrodes 
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Here, C, R, and Jo denote the junction capacitance, re- 
sistance and maximum Josephson current, respectively, 
dots indicate time-derivatives, $o is the magnetic flux 
quantum, and / = Idc + lac sin(a;i) is the total current 
applied to the junction, consisting of a dc and a high 
frequency ac component. The first term on the right 
hand side of Eq. ([1]) describes the displacement current 
CU, where U is the voltage across the junction, and has 
been rewritten in terms of S using the Josephson rela- 
tion S = 27rf//$o- The second term describes the cur- 
rent through the resistor R, the third term the Joseph- 
son current, and the last term the noise current arising 
from Nyquist noise in the resistor. Its spectral power 
density is assumed to be white with Si{f) = iksT/R, 
where T is the temperature and ks Boltzmann's con- 
stant. The model ^ implicitly assumes that magnetic 
fields created by circulating supercurrents can be ne- 
glected (short junction limit). This holds when the lat- 
eral junction dimensions arc below about 4Aj, where 
Aj ~ ($o/47r/xojoAi)^/^ is the Josephson length in terms 
of the critical current density jo, the London penetration 
depth Al, and the magnetic permeability fiQ. 

By integrating Eq. ([T]) one obtains 5 and, by time aver- 
aging, the dc voltage V = <i>o(5)/27r across the junction. 
This is the main observable of our present work, which 
is measured when recording V{Idc), the current voltage 
characteristics (IVC). For numerical simulations, (P) can 
be rewritten in dimensionless units by normalizing cur- 
rents to /o, voltages to IqR, times to tc = $o/(27r/oi?), 
and hence frequencies to fc = IqR/^o, yielding 

/3c'S + 6 + suiS + iN , (2) 

where i = idc + iacSin{Tf / fc) is the normalized applied 
current, r = t/tc the normalized time, f3c = (fc/ fpi)"^ = 
2ttIoR^C /^o the Stewart-McCumber parameter, fpi = 
(/o/(27r$oC))"^^^ the Josephson plasma frequency, and 
iTV the normalized noise current with spectral density 
Si{f / fc) = 4r and noise parameter F = 27rfcBT//o$o- 

In a nutshell, the basic ingredients of NAR as predicted 
in are as follows. The unperturbed deterministic dy- 
namics (Eq. ^ with idc = and i^ = 0) exhibits two 
symmetric attractors, carrying currents of opposite signs 
(zero crossing Shapiro steps). When an external pertur- 
bation in the form of a static bias idc is applied, a subtle 
interplay of this bias force and the dissipation leads to 
a destabilization of that attractor, whose current points 
into the same direction as the applied bias. Its remnant 
is a strange repeller, exhibiting transient chaos, hence 
the name "transient chaos induced NAR" coined in p^ . 
The actual realization of NAR along these lines requires 
a careful choice of model parameters in ([2|) within the 
general regime of frequencies / comparable to fpi and 
values of (3c roughly between 1 and 100. To obtain precise 
quantitative results, we have solved Eq. ^ numerically 
for various such parameter values by integrating and av- 
eraging over typically 5 • 10'^ periods of the ac current. 

For experiments, which were performed at T = 4.2 K, 
we used circular Nb-Al/A10a;-Nb Josephson junctions 



with an area of 200 /im^, cf. upper left inset in Fig. 
1(a). The junctions were shunted by a AuPd strip with 
resistance R = 1.27 Q and integrated in a coplanar wave- 
guide. We denote the critical current Ic as the maximum 
dc current for which = 0. In general, Ic is a function 
of lac and fluctuations. By measuring Ic{Iac ~ 0) and 
matching it with simulations we determined Iq ~ 197 fiA, 
yielding IqR = 250 ^V, fc = 121 GHz and F = 9 • IQ-^. 
The Josephson length is about 40 /xm, i. e., well above 
the 16 fim diameter of our junctions assuring the short 
junction limit. The design value of the capacitance was 
8.24 pF, yielding fpi = 43 GHz, and Pc = 7.9. The ac- 
tual value used in the simulations shown below is some- 
what smaller, namely /3c = 7.7, reproducing particularly 
well the hysteretic IVC in the absence of microwaves. 
The transport measurements have been performed with 
a standard four terminal method, using filtered leads. 
Microwaves between 8 and 35 GHz, with variable output 
power Pm, were applied through a semirigid cable that 
was capacitively coupled to the 50 fl coplanar waveguide. 
The samples were electromagnetically shielded and sur- 
rounded by a cryoperm shield, to reduce static magnetic 
fields. 

Given Iq, R, C, T and Idc, all relevant model param- 
eters are fixed, with the exception of the (frequency de- 
pendent) coupling factor between the microwave ampli- 
tude y^Pm applied from the source and the amplitude 
lac of the ac current induced across the junction. We 
have fixed this factor by comparing the measured depen- 
dence of Ic{\/Pm) with the calculated curve Ic{Iac), as 
shown in the right inset of Fig. [1] for a microwave fre- 
quency of 19GHz (///c ~ 0.16). The experimental and 
theoretical curves are in good agreement. In particular, 
the main side maxima can be found, both in experiment 
and simulation. By adjusting the position of these max- 
ima, we obtain a coupling factor Iac/VPm{iQ GHz) — 
l.OmA/VmW. 

Figure [T]( a) shows IVCs under / = 19 GHz microwave 
irradiation at three values of lac- In the absence of 
microwaves (black line) the IVC is hysteretic, exhibit- 
ing a critical current of 195 fiA and a return current of 
100 A (black arrows). When the microwave field is ap- 
plied, the hysteresis decreases with increasing lac, and 
step- like features appear on the IVC. At Pm = 194 /^W 
{lac = 435 ^A; magenta line), we observe NAR with a re- 
sistance of — 1.07 ri, occurring in an interval \Idc\ < 20 iiA 
(i. e., approximately 10% of /q). When lac is increased 
to P„ = 253 /iW {lac = 497 AtA; green line) the NAR has 
disappeared. However, centered on a voltage which cor- 
responds to the first Shapiro step (Vi = $0/ ~ 39 /iV), 
regions of negative differential resistance appear. In 
Fig. [Hb) measured and simulated IVCs for the two mi- 
crowave amplitudes 435 ^A and 497 fj,A are compared. 
For the former case, which is recorded at the microwave 
amplitude where the maximum NAR has been observed, 
the agreement between the experimental and the theo- 
retical curve is nearly perfect. For the latter case some 
small differences can be seen, although the agreement is 
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FIG. 1; Current voltage characteristics (IVC) of the Joseph- 
son junction at 4.2 K in a 19 GHz microwave field. (a); 
at 3 levels of applied microwave power (0 fiW, 194 fiW and 
253 /iW) showing the effect of negative absolute resistance 
at 194 /iW {lac = 435 nA) and of negative differential resis- 
tance at 253 /iW {lac ~ 497/xA). Left inset: image of the 
Josephson junction. Right inset: measured (thick green) and 
calculated (thin black) dependence of the critical current Ic 
on the microwave current amplitude lac- (b): enlargement 
of the measured IVCs for 194 fiW and 253 fiW, together with 
the simulated IVCs, cf. legend. 

still very good. To demonstrate the origin of the NAR, 
the grey curve in Fig. 1(b) shows a simulated IVC for 
iac = 2.242^^ and F 0, i. c., for the noise-free case. 
The curve shows n = — 1 Shapiro steps to be the cause 
of NAR, clearly revealing its nature to be of the type 
discussed in jlSj . 

Figure [2] compares in more detail the measured 
and calculated dependence of V on Idc and on lac, 
for two frequencies (8 GHz and 19 GHz). For / ^ 
8 GHz, the comparison between measured IcWPm) 
and simulated Ic{Iac) curves yields a coupling factor 
lac/VK'ii&GIiz) = 0.33mA/\/mW. In the graphs, V 
is normalized to $o/, yielding an integer value n for the 
71-th Shapiro step. Again, the agreement between theory 
and experiment is very good. There are at least five lac 
intervals where NAR appears at / = 8 GHz, and three 
such intervals a,t f ~ 19 GHz. Within those regimes, the 
resistance at Idc = reaches values up to about — 1 f2. In 
the case of / = 8 GHz, the NAR persists up to values of 
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FIG. 2: Contour plot of the normalized dc voltage V/^of 
across the junction as a function of dc current Idc and mi- 
crowave current amplitude lac- (a) / — 8 GHz, experiment; 
(b) / = 8 GHz, simulation; (c) / = 19 GHz, experiment; (d) 
/ = 19 GHz, simulation. For symmetry reasons, Idc > —Idc 
implies V i— > —V, hence negative Idc values are not shown. 
Blue areas indicate NAR. 



\Id.c\ ~ 10 /iA, for all values of lac for which NAR shows 
up. In contrast, for / = 19 GHz, the Idc mtcrval for 
NAR decreases with increasing lac- When we increased 
the frequency further to 35 GHz, hysteretic Shapiro steps 
appeared on the IVC, crossing the voltage axis {Idc = 0). 
As a consequence, NAR ceases to exist both in the ex- 
periment and the simulations. 

In a second series of experiments we applied a mag- 
netic field B parallel to the junction plane in order to 
tune (decrease) its Josephson current Iq, making it a i? 
dependent function Io{B) [2^. Thus all /o-dependent 
parameters entering the normalized equation ([2]) aquire 
a _B-dependence, in particular i, Pc, f / fc, and F. Figure 
[3] shows a comparison of the measured and calculated de- 
pendence of the resistance upon lac and Iq{B). Again, we 
find excellent agreement between measurement and the- 
ory. Blue regions indicate NAR. Their most remarkable 
feature is that the values of lac , for which NAR appears, 
practically do not depend on Iq- Furthermore, we find 
that the NAR value can be tuned by Iq via an applied 
magnetic field. For our junction parameters we find a 
maximum NAR at Iq w (0.4 . . . 0.6)/o(B = 0), which is 
increasing with lac- 

In conclusion, we have observed negative absolute 
resistance (NAR) of up to about — Ifi in a shunted 
Nb-Al/A10j;-Nb Josephson junction device subjected to 
microwaves. To clearly see the effect, a careful choice 
of parameters is required, but still the range of suitable 
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FIG. 3: (a) Contour plot of the experimentally measured re- 
sistance R := V{I = 5 fj,A)/5 as a function of the Joseph- 
son current Io{B) and of the microwave amplitude lac- lo 
has been varied by applying a magnetic field to the junction. 
Graph (b) shows the corresponding simulated plot. For both 
graphs parameters at B = are the same as in Fig. 2(a), (b). 



model dynamics reveals that the relevant physical mech- 
anism is of the transient chaos induced NAR type from 
3 ■ The similarity between our Fig. [1] and Fig. 2 in 
9| suggests that with respect to NAR, purely quantum 
mechanical band structure and energy quantization 
effects may be imitated by inertia effects in a purely 
classical, one dimensional noisy dynamics. Moreover, 
our Fig. 3 exhibits many features which are quite similar 
to the corresponding plots in p^ . while the intuitive 
explanation of the almost vertical stripe-pattern in Fig. 
3 remains as an open problem. As an application, our 
present work opens the intriguing perspective of a new 
resistor-type electronic element which is tunable between 
positive and negative resistance via an easily accessible 
external control parameter, e. g., the amplitude of an 
ac driving or an externally applied magnetic field in the 
mT range. 



parameters is quite large. In all cases, we obtain very 
good agreement with theoretical simulations of the 
resistively and capacitively shunted junction model. 
Furthermore a closer inspection of the corresponding 
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